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Abstract. NGC 4815 is a distant and populous open clus- 
ter, which lies in the galactic plane in a region of strong 
absorption. As a consequence, its membership, spatial dis- 
tribution and luminosity function are not well determined. 
In this paper, we present an algorithm which uses both po- 
sitional and photometric data simultaneously to search for 
open cluster members. The contribution of the field stars is 
estimated by our Galaxy model software. The method has 
been developed for deep CCD photometric observations, 
but it can also be used to investigate the membership in 
an n-dimensional space, including proper motions, radial 
velocity and metallicity. 

A clean and well defined colour-magnitude diagram 
(CMD) for NGC 4815 has been obtained after decontam- 
inating from field stars. We have investigated the spatial 
distribution of cluster members and the luminosity func- 
tion. We have shown that the projected surface distribu- 
tion of stars is well-represented by an exponentially de- 
creasing function with a scale length of 1.67 ± 0.06 ar- 
cmin. We have found evidence for mass segregation. The 
luminosity function for NGC 4815 is similar to that of the 
Hyades cluster, and shows a gap at V — 15 mag. 

Key words: open clusters and associ- 
ations:individual:NGC 4815; - Stars: luminosity function, 
mass function - Stars: fundamental parameters - Galaxy: 
stellar content - Methods: statistical 



1. Introduction 

Galactic open clusters can be used to test both the theo- 
ries of stellar evolution (Carraro et al. 1994) and the for- 
mation and kinematics of the Galactic disk (Palous et al. 
1977; Sandage 1988). The advent of CCDs, with high sen- 
sitivity and linearity, has allowed observers to obtain deep 
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photometry of open clusters. This is particularly impor- 
tant for distant clusters or clusters located at low galac- 
tic latitudes and highly contaminated by foreground stars 
(Carraro & Patat 1995). This is the case for NGC 4815 (l 
= 303°, b—-2°), a distant open cluster lying in the galactic 
plane, in a region of strong absorption. 

The first photometric observations of NGC 4815 were 
carried out by Moffat & Vogt (1973), who measured only 
9 stars and did not recognize any sequence. Kjeldsen & 
Frandsen (1991) have obtained UBV - CCD photome- 
try of NGC 4815, including 599 stars, but, because of the 
contamination of noncluster members, they could not de- 
termine cluster membership, reddening or age. More re- 
cently, Carraro & Ortolani (1994, hereafter C094) have 
obtained deeper BV - CCD photometry with the NTT 
3.5 m telescope at the European Southern Observatory 
(ESO). Their sample includes 2498 stars. 

Photometric observations have also been published by 
Phelps et al. (1994) in the context of a survey of old open 
clusters. Unfortunately their observations were made with 
a small telescope (0.9 m) and the errors in the photomet- 
ric data do not allow an accurate determination of the 
physical properties of the cluster. To our knowledge, no 
proper motion study has been published in the literature. 
So, the studies of Kjeldsen & Frandsen (1991) and C094 
will be adopted as the reference photometric studies for 
this paper. 

In order to investigate the cluster properties of NGC 
4815, C094 considered all stars within 4'. 6 of the center 
to be open cluster members. The CMD of their cluster 
members (see Fig. 4 of C094) is very similar to that of 
the total sample (see Fig. 2 of C094). The red stars {B-V 
> 1.5 and V >16) in their CMD of cluster members are 
in fact non-cluster members from the stellar evolutionary 
point of view. Therefore, they cannot reliably discuss the 
spatial distribution and luminosity function of NGC 4815 
because of the many field stars in the sample, and their 
conclusions must be considered as preliminary. 
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Large-format CCD observations can provide deep star 

counts and extend to the edge of the chister, but the larger 
field and fainter magnitude limit also means that cluster 
sequences are contaminated with many field stars. For old 
open clusters, it is well-known that the brighter stars are 
more centrally concentrated (Reid, 1992). Thus, the radial 
distance from the cluster centre is not a suitable criterion 
by which to classify members: by choosing a large radius, 
the contamination of field stars increases, but by using 
a small radius, one can undersample the low mass stars 
which arc not centrally concentrated, and thus bias the lu- 
minosity function. Therefore the development of an algo- 
rithm, which uses all information available (positions, pho- 
tometry, and/or proper motions) simultaneously to search 
for open cluster members, is a necessary task prior to any 
study of many problems in open clusters. 

In this paper, we present a new algorithm from pat- 
tern recognition theory. The advantage of our algorithm 
is twofold. First, we use all available information: appar- 
ent magnitude (V), colour {B ~ V), and positions {X,Y) 
to extract open cluster members in a four-dimensional 
space. Second, we separate open cluster members from 
field stars without assuming any a priori distribution of 
the open cluster members. This is important for our stud- 
ies about the spatial distribution and luminosity function 
of the cluster. We should point out that the method de- 
scribed in this paper is well suited to identifying members 
in a statistical sense, and hence allows one to study the 
statistical properties of the cluster. Additional methods 
clearly need to be used if an analysis requires only cluster 
members. 

The paper is organized as follows: in Sect. 2, we briefly 
describe the structure parameters in our Galaxy model. 
In Sect. 3, we develop a new cluster member finding al- 
gorithm, and investigate the membership of open cluster 
NGC 4815. We have compared our results with those of 
C094 and Kjeldsen & Prandsen (1991). In Sect. 4, the 
spatial distribution of cluster members and the luminos- 
ity function are discussed. Finally, we summarize the main 
conclusions in Sect. 5. 

2. The Galaxy Model 

The difficulty in cluster membership identification is the 
contamination by the field stars. In this paper, we use our 
Galaxy model to predict the distribution of the field stars 
in the direction of the cluster, and establish a classification 
system for identifying cluster members. 

2.1. The starcounts Model 

The main characteristics of stellar populations have been 
extracted from many observations, and the application of 
computer modeling to observations (Bahcall & Soneira 
1980; Pritchet 1983; Gilmore 1984; Robin & Creze 1986; 
Ratnatunga et al. 1989; Reid & Majewski 1993; Chen 



1997a,b) has generated considerable interest. Models can 
parameterize all stellar populations and allow us to calcu- 
late the expected star counts explicitly. 

We have constructed a Galactic structure and kine- 
matic model (Chen 1997a, b), which includes a thin disk, 
a thick disk, and a halo. The model can predict the magni- 
tudes, positions, colours, proper motions, radial velocities, 
metallicities according to the selection criteria used in the 
observation. A detailed description, including kinematical 
parameters (proper motions, radial velocity, space veloc- 
ity), can be foimd in Chen (1997b). In the following lines 
we give a short description of the main structure param- 
eters. 

The method used in our Galaxy model is the numerical 
integration of the fundamental equation of stellar statis- 
tics, for the j-th subsystem; 

A^{V, B-V)=wj ^j{M)Dj{rydr (1) 

Where Aj{V, B — V) is the number of stars of a given 
colour index at a given apparent magnitude; '^j (M) is the 
luminosity function and Dj{r), the density law. In order 
to predict star counts in a given direction, a luminosity 
function and density distribution for each population or 
sub-population in the model is required. The stellar den- 
sity laws used for the thin and thick disk are exponential. 
For the population II spheroid stars, the density law is 
adopted following de Vaucouleurs (1977). We adopt the 
luminosity function for the thin and thick disk stars given 
by Wielen et al. (1983). The halo luminosity function is as- 
sumed to be similar in shape to that of the disk for My > 4, 
but matching a globular cluster function at brighter mag- 
nitudes (Reid, 1993a). We adopt a CMD given by Bahcall 
et al. (1987) from the Yale Parallax catalog for thin and 
thick disk stars. For halo stars, the colour-magnitude di- 
agram for both the main-sequence and evolved (red giant 
and horizontal-branch) stars is based on the table given 
by Bergbush & Vandenberg (1992). 

The thin and thick disk giants were included separately 
in the model. According to Bahcall & Soneira (1981), we 
take the fraction of stars, /, on the main sequence (MS) 
in the plane of the disk to be: 

/ = o.44e"-OOOi5(M„+8)^-^^ My < 3.7 

f = l,My> 3.7 

We have allowed for observational uncertainties by 
adding Gaussian errors to the V magnitude and to the 
colour of each 'star' generated. C094 have obtained er- 
rors of 0.02, 0.05, and 0.07 mag at ^=15.5, 17.5 and 19.5 
mag respectively, and at the same magnitude levels, they 
found dispersions in colour of 0.05, 0.11 and 0.14 mag, re- 
spectively. These results were incorporated in our Galaxy 
model software. 

We should point out that these parameters in the 
model code are not the results of this investigation, but 
have been derived from previously published results (Bah- 
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call & Soneira 1980; Robin & Creze 1986; Reid & Majew- 
ski 1993; Chen 1997a; Chen 1997b). 



2.2. Extinction in the line of sight 

One of the important characteristics of galactic obscura- 
tion is that it is patchy, so, it is important to discuss the 
foreground obscuration in the region of NGC 4815 {d ^ 
2300 pc). Neckel & Klare (1980) have derived extinctions 
and distances for more than 11000 stars and investigated 
the spatial distribution of the interstellar extinction at | h 

1 < 7°. 6. Their results are shown in a series of diagrams. In 
Fig. 1, we plot their results (open circles) from the nearest 
field (/ = 302°, h = -3°) in the direction of NGC 4815. Us- 
ing a large sample (about 17000 stars) with MK spectral 
types and photoelectric photometry, Arenou et al. (1992) 
constructed a tridimensional model of the galactic inter- 
stellar extinction. In Fig. 1, we show the results (solid line) 
derived by Arenou et al. (1992). 

Recently, we have used Lund Catalogue of Open Clus- 
ter Data (Lynga 1987) to derive an analytic expression of 
the interstellar extinction (Chen et al. 1997b). In Fig. 1, 
we have overplotted the average extinction in the galac- 
tic plane (filled circles). We can see that extinction laws 
derived from different methods are fairly consistent. More- 
over, Bertelli et al. (1995) derived the extinction in the line 
of sight in the field near the cluster Lynga 7 {I = 328°, b 
= -2°. 8), not far from the field of our investigation. They 
found an extinction of Ay = 1.5 mag at 2000 pc, and Ay = 

2 mag at 4000 pc, which is also consistent with the results 
shown in Fig. 1. 

In Fig. 2, we show the observed CMD for all the stars 
in a field of view of 12'.3 x 12'.3 from C094. We can 
recognize the MS and the giant branch. Many field stars 
seem to contaminate the CMD, especially a red sequence 
parallel to the MS. 

In Fig. 3, we show the simulated colour-magnitude dia- 
gram from the Galaxy model with the Arenou et al. (1992) 
extinction law. This predicted diagram provides us with 
the distribution of field stars in this direction. In following 
discussion, we adopt the extinction law derived by Arenou 
et al. (1992) in our Galaxy model. 

3. Identification of NGC 4815 members 

Many authors have used colour-magnitude and colour- 
colour diagrams to separate open cluster members and 
field stars. When kinematic data are available, it is com- 
monly accepted that the membership probabilities ob- 
tained from the analysis of proper motions or radial veloci- 
ties are more reliable. The parametric models proposed by 
Vasilevskis et al. (1958) and Sanders (1971) assume that 
the proper motion distributions of field stars and cluster 
members can be modeled by a normal bivariate function. 
However, in several cases the wings of the field star dis- 
tribution are more extensively populated than expected 
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Fig. 1. The extinction along the line of sight in NGC 4815. 
Solid line is the Arenou et al (1992) extinction law, filled circles 
arc the results derived by us (Chen ct al. 1997b) from open 
cluster system, open circles are the results of Neckel & Klare 
(1980). NGC 4815 is located at about 2300 pc from the Sun. 
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Fig. 2. Observed CMD for all the stars studied in the region 
of NGC 4815 
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Fig. 3. Simulated CMD in the region of NGC 4815 



based on the predictions of a normal bivariate function 
(Marschall & van Altena, 1987). Cabrera-Cano & Alfaro 
(1990) have developed a non-parametric approach to iden- 
tify open cluster members by using proper motions and 
positions as variables. However, because their method as- 
sumes that the variables are independent, it cannot be 
used with photometric data. 

In this paper, we have developed a new cluster mem- 
ber finding algorithm. This method, based on a non- 
parametric approach, has been previously used to identify 
moving group members in the solar neighborhood (Chen 
et al. 1997a) and to isolate galactic thick disk stars from 
star count survey (Chen et al. 1992; Chen 1996a). The vi- 
ability of the method has been tested by a series of Monte 
Carlo simulations. Results show that the method is pow- 
erful. The detailed descriptions and numerical simulations 
of the method can be found in Chen (1996b) and Chen et 
al. (1997a). 



3.1. The method 

Let (xi, X2, X:i, ... Xn) be a set of stars to be studied. Each 
star Xi is considered as one point in the four-dimensional 
space {X, Y,V,B — V) of position, magnitude and colour. 

Let Ni and be the number of stars in the observed 
and the 'simulated' sample, respectively. The number of 
stars {N2) comes from the Galaxy model and is considered 
as the contribution to the counts of the field stars in the 
direction of the open cluster. According to the Bayesian 



theory, the probability of a given star being an open clus- 
ter member can be derived as: 



P{c\Xi) = 1 - P{f\Xi 



N2p{x,\f) 
Nip{x,) 



(2) 



Where P{c\xi) and P(f\xi) are the a posteriori prob- 
ability that an object Xi belongs to the cluster and field 
stars, respectively. p{xi) is the probability density func- 
tion (hereafter pdf) from the observed sample, p{xi\f) is 
the class conditional pdf for field stars. The probability 
P{c\xi) can be derived by calculating p(a;i) a,nd p{xi\i) at 
each observed star Xi. 

The true p{xi) in each star can be estimated by a ker- 
nel estimator p (xi) of the density (Hand, 1982): 



1 



|E|'/'(27r)<i/2 



(3) 



where Xi is the point at which the estimate is being made, 
n is the number of stars in the observed sample, xj is the 
observed sample set defined in d-dimensional space, {xi 
-Xj)' is the transpose of the vector (xi -Xj) and E is the 
variancc-covariance matrix of the observed sample. 

The optimal smoothing parameter h derived by Silver- 
man (1986) can be written as: 



'd + 2' 



l/(d+4)^„-l/(d-F4) 



(4) 



where a is the average marginal variance, — d^^Scrf . 

From the definition of the kernel estimator, we know 
that it is a sum of the multivariate normal density func- 
tions placed at the observations. The functions determine 
the shape while the smoothing parameter h determines 
their width. Concerning the computation of the class con- 
ditional probability density function, p{xi \ /), the proce- 
dure is the same as described above for p{xi), but we have 
to perform the summation over the simulated sample for 
the field stars. 

3.2. Separation of field stars from cluster members 

The method described above has been used to study the 
membership of the open cluster NGC 4815 in the four- 
dimensional space of {X, Y,V, B — V). C094 investigated 
the degree of completeness for their observations, and 
found that the sample is almost complete to F = 17 mag. 
In order to avoid incompleteness problems, we only ana- 
lyze the stars with this limit, thus including 739 stars. The 
determination of those stars considered as cluster mem- 
bers has a probabilistic nature. In Fig. 4, we show the 
CMD for cluster members with probability P{c\xi) > 0.5, 
0.8, and 0.9, respectively. Isochrones (Alongi et al. 1993) 
of different metal abundance Z and age are overplotted. 

In order to consider the influence of the observational 
errors {ay and as-v) in our analysis, we have gener- 
ated several sets of 'observed' data by adding an incre- 
ment {AV, and A{B — V)) to the real values of each star. 
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Fig. 4. Left panel: The CMD for the cluster members, the Arenou et al. (1992) extinction law has been adopted. Superimposed 
are isochrones of different metal abundance Z. Dashed line represents a Z= 0.008 isochrone for an age of 500 Myr, plotted 
adopting Eb-v = 0.70 and {V — Mv) = 14.2. Solid line represents a, Z = 0.004 isochrone for an age of 600 Myr, plotted adopting 
Eb-v = 0.75 and {V — Mv) — 14.10, dotted line represents a, Z — 0.02 isochrone for an age of 500 Myr, plotted adopting Eb-v 
= 0.65 and {V — My) = 14.20. a) The cluster members with probability P{c\xi) > 0.5, b) The cluster members with probability 
P{c\xi) > 0.8, c) The cluster members with probability P{c\xi) > 0.9 

Right panel: The CMD for the field stars, d) The stars with probability P{c\xi) < 0.5, e) The stars with probability P{c\xi) < 
0.8, f) The stars with probability P{c\xi) < 0.9 



Each of these increments has been randomly chosen from 

a normal distribution of zero mean and standard deviation 
equal to the individual observational error. We found that 
the member probabilities are not significantly changed, 
indicating that the observational errors are not very im- 
portant in the present analysis. 

As mentioned previously, we have used the Arenou et 
al. (1992) extinction law in our Galaxy model. In order 
to consider the influence of the uncertainty due to the ab- 
sorption in our analysis, we have checked our results by 
using the Arp (1965) reddening distribution: the absorp- 
tion increases at about a rate of 1.0 mag kpc~^ for the first 



2 kpc, and after 2 kpc, there is no further light absorp- 
tion. This simple reddening law has been found to be in 
agreement with the observations along several galactic di- 
rections in the galactic plane (Arp, 1965; Paczynski et al. 
(1994), Ng & BertelH, 1996; Arp & CufFey 1962). In Fig. 5, 
we show the results from the Arp reddening distribution. 
In Fig. 6. we show the P(c|a;j) from the Arenou model (PI) 
against the P{c\xi) from the Arp reddening distribution 
(P2). ^From Fig. 6, we can see that, for large PI and P2 
(> 0.8), the difference between PI and P2 is small, but 
PI is systematically larger than P2; for small PI and P2 
(< 0.5), the difference between PI and P2 becomes larger. 
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Fig. 5. as Fig 4, but using the Arp (1965) reddening distribution. 



For example, with P{c\xi) > 80%, we found 378 stars with 
the Arenou et al. (1992) extinction model, and 314 stars 
with the Arp reddening distribution. Therefore, we can say 
that the uncertainty due to the absorption in our analysis 
does not change too much the likely open cluster members 
(P(c|a;j) > 80%). However, we should point out that, to 
derive more reliable membership probabilities, one should 
have a good knowledge about the extinction law in the 
direction of the cluster. 

To our knowledge, the only investigation in the liter- 
ature of the cluster membership for NGC 4815 is that of 
C094. They separate cluster stars from field stars by us- 
ing a radius of 4'. 6 arcmin. In Fig. 7, we show the CMD 
for the cluster members by the criterion of C094 for stars 
with V < 17 mag. It can be seen that the use of only 
position to determine the cluster membership leads to a 
significant contamination by field stars {V > 16 and S — V 
> 1.5). In a standard cluster analysis, identifying cluster 



members is a process of attrition - first selecting by posi- 
tion (or proper motion), then removing photometric non- 
members, then spectroscopic non-members and so forth. 
Our method, combining position and (V. B — V) into a sin- 
gle step, enlarges the statistical "distance" between the 
field stars and cluster members, thus reducing the con- 
tamination due to the field. Moreover, our method, using 
a probability classification, can provide an unbiased esti- 
mation of the intrinsic distribution of the cluster stars. 

Another database for NGC 4815 in the literature is 
that from Kjeldsen & Frandsen (1991), who obtained 
UBV CCD observations. In Fig. 8, we show the colour- 
magnitude diagram for their observations. As pointed out 
by Kjeldsen & Frandsen (1991), we can see that the CM 
diagram contains many non-members. In any case a MS 
with a turn-off is still clearly visible. We found that the 
isochrone of Z = 0.008 for an age of 500 Myr can basically 
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Fig. 6. The probabihty PI derived from the Arenou et al. 
(1992) extinction model against the probability P2 derived 
from the Arp (1965) extinction model 



Fig. 8. CMD obtained by Kjeldsen & Prandsen (1991). Super- 
imposed is isochrone oi Z = 0.008 for age — 500 Myr, plotted 
adopting Eb-v= 0.7 and V - M„=14.2 mag 
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Fig. 7. CMD of NGC 4815 for cluster members derived from 
the criteria of C094, i.e. stars within a radius of 4.6 arcmin of 
cluster center. 



Fig. 9. The luminosity function of the cluster NGC 4815. The 
luminosity functions of the Hyades (filled circles) and NGC 
7789 (open circles) , which are shifted arbitrarily to match the 
bright part, have been overplotted 
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Fig. 10. The Mass function of the cluster NGC 4815. The sohd normalized to the surface density at r= 1.2 arcmin. 

line is the result from least-squares fit 

explain the observations of Kjeldsen & Frandsen (1991), 
which is consistent with our result from C094 sample. 
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Fig. 11. The radial density distribution for the cluster NGC 
4815, where the error bars were calculated assuming that the 

number of stars in a bin is governed by Poisson statistics, the 
dotted line we have adopted to fit the spatial distribution is 
an exponentially decreasing function of the distance from the 
cluster center with a scale length of 1.67 arcmin 



4. The spatial distribution and luminosity function 

of NGC 4815 

In this section, the cluster member probability discussed 
above is used to investigate the spatial distribution and 
luminosity function of NGC 4815. Both C094 and this 
investigation have shown that NGC 4815 has an age of 
about 500 - 600 Myr, which is similar to the age of the 
Hyadcs cluster (625 ± 50 Myr) derived by Lebreton et al. 
(1997) from Hipparcos observations. It is very interesting 
to compare our results with that of the Hyades cluster, 
which has been widely studied and whose properties are 
well understood. 

4-1- Luminosity function 

It is very important to study the cluster luminosity func- 
tion because it can provide information about both the 
initial mass function and the cluster dynamical evolution. 
Cluster membership probabilities were used to determine 
the luminosity functions for the cluster. The luminosity 
function $(Mt,) was obtained as: 



(5) 



Where P{c\xj) is the membership probability of star j, 
AMy is the bin size. Table 1 and Fig. 9 show the resulting 
luminosity function for NGC 4815. 
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Our results show that the luminosity function increases 
aiV = 14.0 mag (M^, = -0.2) around the MS turnoff, with 
a prominent drop at y = 15 mag (M„=0.8) that corre- 
sponds to the MS gap. The presence of a gap in the dis- 
tribution of stars near the MS turnoff region is a common 
feature of old (M67, Montgomery et al 1993) and inter- 
mediate age (NGC 7789, Roger et al. 1994) open clusters. 
The gap is considered as a real feature corresponding to 
the evolutionary behaviour of stars at the end of their 
core H-burning stage when they undergo the overall con- 
traction phase (Carraro et al. 1994). This gap can also be 
seen in the CMD for all observed stars as well (see Fig. 2). 
However the luminosity function derived by C094 did not 
show this feature, probably due to the presence of many 
field stars that were mis-classified as cluster members. In 
Fig. 9 we have overplotted the luminosity function of the 
Hyadcs clusters from Reid (1993b) and NGC 7789 from 
Roger et al. (1994). 



Table 1. The luminosity function of the cluster NGC 4815 
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Given the luminosity function derived from our data, 
we can apply mass-luminosity relations to derive an esti- 
mate of the total mass of the NGC 4815. We have used 

the relation 

log{M) = -QMlMy + 0.356 M„ < 0.22 (6) 

log{M) = -0.163M^ -|- 0.558 My > 0.22 (7) 

derived from Andersen's data (1991). ^Prom masses 
[Ai, in solar mass) and membership probabilities, the ob- 



served mass in the cluster can be determined from T,j 
P(c|a;j)A^I . To the limit of the photometry, wc find a mass 
of the cluster of 912 solar masses. In Fig. 10, we show the 
derived mass function {N{log{M))) for the NGC 4815. By 
using a least-squares fit, we found that the mass function 
has a slope of -2.2 ± 0.3 for 1.5 < < 2.3 in the plane of 
{log {N), log{A4). Francic (1989) has derived the slope for 
the composite mass function for six clusters with a mean 
age of 400 Myr, he found d log N/d log (M) = -1.97 ± 
0.17 for 1.1 < Ai < 2.5. His result is in agreement with 
our result in this paper, considering the error bars. Our 
result is also in good agreement with the value determined 
by Scalo (1986) for the field stars (-2.0 ± 0.18 for 0.8 < 
M < 18). 

4-2. The spatial distribution of cluster members 

We have taken annuli with a width of 0.3 arcmin centered 
on the cluster center and have counted the numbers of 
stars {N = P{c\xj) falling within each annulus. Fig. 11 
shows the radial density distribution for the cluster NGC 
4815, where the error bars were calculated assuming that 
the number of stars in a bin is governed by Poisson statis- 
tics. 

It was shown by Van den Bergh & Sher (1960) that the 
projected surface distribution of stars for an open cluster 
can be represented by p{r) = poexp{—r/ro), where tq is 
the scale length for the projected surface density of clus- 
ter members. We have used this exponentially decreasing 
function of the distance from the cluster center to fit our 
results, and found the scale length to be ro = 1.67 ± 0.06 
arcmin. Since the distance to NGC 4815 is known, we can 
transform the angular scalelength to a linear distance. As- 
suming that the cluster is symmetric, wc can calculate the 
average volume density of stars within the cluster. For the 
inner region of the cluster (r < 1.67 arcmin), we derive an 
average density of 47.5 stars/pc"^ and a mean separation 
of 0.171 pc. If we calculate the density at the edge of the 
cluster, taken to be at a diatance of 4.0 arcmin (see Fig. 
11). we find an average density of 7.9 stars/pc^ with a 
mean separation of 0.312 pc. 

Reid (1992) has studied the Hyades cluster members 
and found that the brighter stars arc more centrally con- 
centrated (mass segregation). In order to look for similar 
evidences in NGC 4815, we have divided our cluster mem- 
bers into three magnitude ranges, 15 > V > 14: 16 > ^ 
> 15; 17 > y > 16. In Fig. 12 we plot the surface density 
distributions for these three magnitude ranges, all distri- 
butions have been normalized to the surface density at r= 
1.2 arcmin. From Fig. 12, we can see that the brighter stars 
are slightly more centrally concentrated, as one would ex- 
pect from mass segregation. However, because the sample 
used here is limited with V < 17 mag, the observational 
evidence is limited. Deeper samples are needed to confirm 
this result. 
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5. Conclusions 

We have developed a new algorithm, which uses both po- 
sitional and photometric data simultaneously to search for 
open cluster members. The algorithm can be extended to 
use all possible variables, including proper motions, radial 
velocity and multi-color photometry automatically. When 
the contamination of field stars is not very serious, then 
simple methods can be used to identify cluster members 
(radial distance from the centre; position in the colour- 
magnitude diagram; etc); if field star contamination is sig- 
nificant, then one requires a more sophisticated technique. 
The method developed in this work, which uses all vari- 
ables simultaneously to search for open cluster members, 
gives better results. This is particularly true when cluster 
stars lie beyond the observed sample that is used, where 
the contamination of field stars is difficult to estimate. 

The algorithm has been used to isolate open cluster 
members of NGC 4815 from the field stars. We have com- 
pared our results with those of C094. A clean and well 
defined colour-magnitude diagram (CMD) for the cluster 
has been obtained from our new algorithm. After deter- 
mining the cluster member probability, for the first time 
we have investigated the spatial distribution of the cluster 
NGC 4815 and its luminosity function. To summarize, we 
find that: 

1. NGC 4815 is a cluster of the Hyades generation, with 
an age of about 500 - 600 Myr. 

2. the projected surface distribution of stars is well- 
represented by an exponentially decreasing function 
with a scale length of 1.67 ± 0.06 arcmin. 

3. there is evidence for mass segregation. 

4. the luminosity function for NGC 4815 is similar to that 
of Hyades cluster, and shows a gap at F = 15 mag. 

5. to the limit of the photometry, the observed mass is of 
about 900 solar masses for NGC 4815 

6. the mass function has a slope of -2.2 ± 0.3 for 1.5 < 

< 2.3 in the plane of {log (N), log (M)). 

These observational results are important, and can be 
used to constrain models of stellar evolution (Carraro et 
al. 1994). Other contaminated open clusters are going to 
be analyzed in the future. 
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